Thermal entanglement in fully connected spin systems and its RPA description 
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We examine the thermal pairwise entanglement in a symmetric system of n spins fully connected 
through anisotropic XY Z-type couplings embedded in a transverse magnetic field. We consider both 
the exact evaluation together with that obtained with the static path + random phase approximation 
(RPA) and the ensuing mean field + RPA. The latter is shown to provide an accurate analytic 
description of both the parallel and antiparallel thermal concurrence in large systems. We also 
analyze the limit temperature for pairwise entanglement, which is shown to increase for large fields 
and to decrease logarithmically with increasing n. Special finite size effects are as well discussed. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 



I. INTRODUCTION 

Quantum entanglement, one of the most fundamen- 
tal and intriguing features of quantum mechanics, is well 
recognized as an essential resource for quantum informa- 
tion processing and transmission [H-Q. It has recently 
acquired an important role also in many-body and con- 
densed matter physics 0-0| , where it provides a new per- 
spective for analyzing quantum correlations and quan- 
tum phase transitions, as well as in other fields like the 
foundations of statistical mechanics [8( . The study of en- 
tanglement in intcractingspin models has in particular 
attracted much interest since they provide 

a basic scalable qubit representation suitable for imple- 
menting quantum processing tasks and are at the same 
time able to capture the main features of diverse physical 
systems. Some of these models can in addition be exactly 
solved for any size, providing hence a suitable scenario for 
testing the accuracy of approximate descriptions. 

An example is that of a symmetric array of n fully con- 
nected spins (simplex) with anisotropic XYZ type cou- 
plings embedded in a uniform transverse magnetic field. 
This is a solvable yet non-trivial model which exhibits a 
quantum phase transition at T = 0, whose Hamiltonian 
is formally equivalent to that of the well-known Lipkin- 
Meshkov-Glick (LMG) model PH. It has attracted re- 
newed interest in recent years, having been used to de- 
scribe diverse physical systems such as Josephson junc- 
tion arrays [l9[ and two mode Bose Einstein condensates 
[2fj| . Its zero temperature entanglement properties were 
analyzed in detail in refs. [12H1 5j | . where it was shown 
in particular that the pairwise concurrence, a measure 
of the entanglement between two spins [2l[ , exhibited a 
rich behavior when properly scaled, with a cusplike max- 
imum at the critical field in the ferromagnetic case and 
a smooth decrease for large fields [3]. 

In this work we will examine the thermal pairwise en- 
tanglement in this system, together with its description 
in the framework of the mean field + random phase ap- 
proximation (RPA) derived from the path integral rep- 
resentation of the partition function [22l |23[. Our aim 
is twofold. First, we want to determine its thermal be- 
havior and stability, a relevant aspect in physical real- 



izations. Secondly, given the complexity of determining 
the entanglement properties in general interacting many- 
body systems at finite temperature, we want to examine 
the extent to which its main features can be captured 
by a general tractable method like the RPA, which takes 
into account just small amplitude quantum fluctuations 
around the mean field. We will show that the present 
RPA treatment provides, for anisotropic couplings, an 
accurate analytic description of both the parallel and 
antiparallel thermal pairwise entanglement in large sys- 
tems, generalizing the results of [23| for the XXZ case 
(where entanglement is just antiparallel). The limit tem- 
perature will be shown to decrease only logarithmically 
with increasing n at all fields, for the standard \ jn scal- 
ing of coupling strengths, and to exhibit a different field 
dependence in the parallel and antiparallel sectors. In 
particular, it increases for large increasing fields, despite 
the decrease of the T = concurren ce ( and at variance 
with the behavior in the XXZ case _23|), entailing just 
a finite separable field window at any temperature. 

Section II describes the model and its exact partition 
function and concurrence, together with their evaluation 
in the static path and mean field + RPA treatments and 
the asymptotic expressions. Section III discusses in de- 
tail the exact numerical and approximate results in finite 
systems. Finally, conclusions are drawn in IV. 



II. FORMALISM 

A. Exact partition function and concurrence 

We will consider n qubits or spins 1/2 coupled through 
an anisotropic full range XYZ Heisenberg interaction in 
a transverse magnetic field b. The Hamiltonian reads 

i=i ijtj 

= bs z -l £ » M (s2-5), (1) 



H=x,y,z 



where denotes the spin component at site i (in units of 
h) and = y^j—i s 1 ^ the total spin components. The 1/n 



2 



scaling of coupling strengths ensures that all intensive 
energies remain finite for n — > oo and finite w M . For T > 
the total spin S 2 = ^TJ Sft is n0 longer fixed, so that all 
terms in |T]) are independent. Nonetheless, without loss 
of generality we can assume \v y \ < \v x \ and b > 0. We will 
consider here the attractive case v x > (with \v y \ < v x ) 
where the ground state will have maximum spin S = n/2. 

Since H is completely symmetric and commutes with 
both S 2 and the SVparity P — exp[i7r(S' z + n/2)] (global 
phaseflip) the partition function at temperature T = 
(we set Boltzmann constant k = 1) can be written as 



s/2 



Z = Tiexp[-(3H}= J2 Y ( S ) Yl e- f3Esk 



(2) 



S=5 n 



u=±,k 



where Y(S) = ( n/ »_ s ) - ( n/2 » with Y(n/2) = 1, 

is the multiplicity of states with total spin S, such that 
J2s=s n Y(S)(2S + 1) = 2™ [S n = (|) for n even (odd)] 
and Esk" are the eigenvalues of H with total spin S and 
parity v [k = 1, . . . , ff+g+Wg — #»,)]■ It should be noticed 
that in the fermionic realization [181 ] . the multiplicities 
Y(S) would be different (the total number of states in 
the half-filled fermionic system is ( 2 ™) instead of 2"). 

The pairwise entanglement at T > is determined 
by the reduced two-spin density matrix pij = Tr n^^yp 
(i 7^ j), where p = Z^ 1 cxp[— j3H] is the global thermal 
density, will be entangled if it cannot be written as a 
convex combination of product densities pij ]. i.e., if pij 7^ 
J2 a 1aPi <8> p" , with g Q > 0, being separable otherwise. 
The amount of pairwise entanglement can be measured 
through the entanglement of formation which 
in the case of two qubits can be evaluated as [2l[ = 

~ Y,u=± 1v lo §2 with 1± = K 1 ± ^/l - Cfj) and Cy 
the concurrence [2l| , itself an entanglement measure [26j j . 
Since By- is in this case just an increasing function of 



Cij, with Cjj 



Eij — 1 (0) for a maximally entangled 
(separable) pair, it is equivalent to use Cy as measure. 

In the present system p^ will be the same for any pair 
and will commute with the reduced parity exp[i7r(s* +s° z + 
1)] and total spin X^( s m + s i,) 2 i being in the standard 
basis of the form 



Pij = 







( p+ 

po a_ 

a_ po 

V a+ p_ 



P± : 
Po = 



(44; 
- a z 



4 



where s!(_ = s x ± is^ and (p — x, y, z) 



T <91nZ 



" n-1 dv M 

TdlnZ 



n db 



(i * J) , (3) 



(4) 



Note that - 3 (^ T y < a M < \ as (S 2 ) = f + n(n - 

The ensuing concurrence C = can be expressed as 



C = Max[C + ,C*_,0], with 

C + = 2(|aq_| - po) = 2(a x - a y + a z - 
C-=2(\a_\-^pTp-) 

= 2(a x + a y - ^(L + az )2- {Sz )2) 



(5) 
(G) 



Here C+ (C ) denotes a concurrence of parallel (antipar- 
allel) type [23], as in Bell states |tt) ± \U) (\U) ± lit))- 
Just one can be positive for a given pij. In the final 
expressions (J!)])-® we have assumed ct± > 0, valid for 
the present attractive case v x > \v y \. Since all pairs are 
equally entangled, the maximum value that can be at- 
tained by C in the present system is 2/n [28[ (reached 
for instance in the W-state \SM) — |f , § — 1)) implying 
that only the rescaled concurrence c = nC can remain 
finite in the thermodynamic limit n — > 00. 



B. Static path + RPA 

The auxiliary field path integral representation of the 
partition function ^ can be written as [22] 







/ D[r]Tr 


Texp{- / H[r(T)]dr} 




Jo 



H{r) =bS z -r-S 



1 EC 



(7) 



(8) 



where r = (x,y,z), T denotes (imaginary) time or- 
dering and H(r) represents a linearized Hamiltonian. 



nr (t) 



dr] 



1 is as- 



Normalization / D[r] exp[- f Q ^ M 4lv 
sumed. Starting from a Fourier expansion r(r) = 
EZ-oo r k e iUkT , uj k = 2nk/p, with D[r] cx IL^fc, the 
static path + random phase approximation [291 - 4311 ] [to be 
denoted as correlated SPA (CSPA)] preserves the full in- 
tegral over the static components r = r$ but integrates 
over r*fc, k =/= 0, in the saddle point approximation, for 
each value of the running static variables. This proce- 
dure takes thus into account large amplitude static fluc- 
tuations, relevant in critical regions, together with small 
amplitude quantum fluctuations, and is feasible above a 
low breakdown temperature T*. The final result for the 
present spin 1/2 system can be cast as 

^ A(r) sinh[i/3o;(r)] 
{x,y,z- b), 



^CSPA — 



n(3 
Attvh 



where, defining A 
Z(r) 

A(r) 
u(r) 



b 

Trexp[- (3 H(r)} 

e -iP E M («</^+^) [ 2 C o S h §£A(r)]" , 

iJ2^-u>)(i-U")] 1/2 > 

v„ tanh[i/3A(r)]/A(r), 



(10) 
(11) 

(12) 
(13) 
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with ix' < /j/,fi" 7^ /it. In © Z(r) is a Hartree-like 
partition function while the remaining factor accounts 
for the small amplitude quantum corrections, with uj(r) 
the single collective thermal RPA energy existing in the 
present system. It can be obtained from the equation 



DetfJ,, 



E 



V-V - Pu , 



0, (14) 



with s£ = (i/|s m |-i/), p v = e P e "/Y,u e 0e " and k)) £ ^ 
the eigenstates and eigenvalues of A-s. If u M < 0, the cor- 
responding integral should be done along the imaginary 
axes and can be evaluated in the saddle point approxi- 
mation 3l|. The elements ©-(jU) become 



1 



1 ' coth±/3u; — 



and (s z ) = h{z)/v z , where (. . .) denotes CSPA averages. 



C. Mean field + RPA 

For sufficiently large n and away from the critical re- 
gion, we may integrate all variables r^, including ro, in 
the saddle point approximation around the minimum of 
the free energy potential —ThxZ(r), determined by the 
self-consistent equations 



M> IJ> = x,y,z. 



(15) 



This leads to the mean-field+RPA (MF+RPA). For an 
isolated minimum at r = ro , we obtain 



-'MF+RPA — 



Z(r ) sinh±/3A 
VI — C sinh hj3w 



(16) 



where C, = 1 - ^Det[-^ ^r^§^ }r accounts for the 
gaussian static fluctuations and A = A(r ), w = w(r ). 
In (fT6|) . Z(ro) is the MF partition function while the last 
factor is the proper RPA correction, which represents the 
ratio of two independent boson partition functions: that 
of bosons of energy to to that of bosons of energy A. 

For the present Hamiltonian Eqs. (|15p imply either 
r M = or = 1 for fi = x,y. For \v y \ < v x and v z < v x , 
we then obtain the following minima: 
a) If |6| < b c and T < T c (b), where 



Tc(b) 



m l-b/b c 



(17) 



the minimum corresponds to the degenerate parity- 
breaking solution r = (±x, 0,z), with x =/= 0. In this 
case A is determined by the equation f x = 1, i.e., 

A = v x tanh ±/3A , (18) 

which depends just on v x and T (A = v x at T = 0), while 
z = —v z b/b c (independent of T) and x = \J\ 2 — v 2 b 2 /b 2 , 



the constraint A > v x b/b c leading to Eq. (IT7|) . At this 
solution, the RPA energy (IT21 becomes 



w = s,/(l -/«)(!-/.). 



(19) 



with ffj — v^/v x , while £ = h(3v x / cosh i/3A. Note that 
w ->• for T -> T c (6) (as a; -> 0) or n s -> v x (as / y -» 1), 
implying the divergence of (fTB")) in these limits (see [23[ 
for the correct MF+RPA treatment in the continuously 
degenerate XXZ case). 

b) For |fe| > b c or T > T c (b), the minimum corresponds 
to the normal solution r = (0, 0, z). In this case A = b — z 
is the positive root of the equation 

X = b + v z tanh[|j8A] (20) 

with A = b + w 2 at T = 0. The RPA energy becomes 



V(l -/*)(!-/„)■ 



(21) 



with / M = (1 - b/X)v^/v z , while C = i/3w 2 /cosh 2 ±/3A. 
Here w — > for T — > T c (b) (as f x — ► 1) but remains finite 
for v y — > v x - This is also the only solution for v z > v x . 
The ensuing expressions for the elements ©-(J!]) are 



in ., 



2(n-l) v 2«2 2 



1 ,nz . , 
: 2^ ( - + ^' 



9A ,, -, duj ,. i _ 

^ = "S - coth ~ "S - coth 5P W + T > o i 
or] or) 1 — Q or] 

with ?7 = v^,b. The first term in a^, (s z ) is the O(l) 
Hartree contribution, whereas 8^ provides the 0(1/ n) 
RPA corrections, essential for describing entanglement. 

D. Asymptotic expressions for the concurrence 

Full expressions for the MF+RPA concurrence are 
rather long and are given in the Appendix. However, 
up to 0(1/ n) terms and for sufficiently low T, we obtain 

C+«— ^-(1 — cothi/3o;)-2e-' 3 ^, (22) 



71—1 V x - Vy 



' ' ? ^^cothi/3w)-2 e -^ (23) 



-(1- 



n — 1 ui 
in the symmetry-breaking phase (|6| < b c ), where 



'l-(b/b c y 



X 



Vy - V z 



(24) 



JX Vy 



' i-x ' 

whereas in the normal phase (b > b c ), C- < while 

C+ = -^—(1 cothi/3w) - 2e-^ b+v $5) 

n-V b + v z -v y 2 ' y ' 

with 



b/b c - 1 



b + v z -v y V b/b c - x 



(26) 



4 



We have used in (|22 ]) -([26 |) the T = values for A and 
lj, as the ensuing thermal corrections will be of order 
e~^ A = 0(l/n) for temperatures where |22|) . (f23|) and 
([23)) are positive, leading then to 0(l/n 2 ) terms in C. 
Eqs. ([22]) (|25| become increasingly accurate as n in- 
creases (coinciding for T — > and w z = with the ex- 
pressions of ref. [14f ) and can be summarized as 



l-^^cothi^ 



i 



2e 



-/3A 



(27) 



with A = w = -^(1 - 6 2 /6§)(A — ^JCA — for \b\ < 
b c and A = b + v z , w = — v x )(X — v y ) for b >b c , the 
result for b > b c applying just for C+. 

For T -> (e- /3A -> 0, coth ±/3w -> 1) C± is then fully 
determined for large n by the scaled field b/b c and the 
anisotropy \- F° r < X < 1 (i- e -j v z < v y < v x ), C_ 
(C+) will be positive for |6| < b s (> b s ), where 

b s =b c Vx, (28) 

is the factorizing field [2 .32] , where the system possesses 
a separable ground state. Accordingly, at T — both 
C± vanish at b = b s , C being antiparallel for \b\ < b s 
and parallel for |6| > b s . On the other hand, if x < 
(v y < v z < v x ) or x > 1 («» < < v z , in which 
case there is no symmetry-breaking phase) C is always 
parallel at T = 0. It is also seen from (|25p - (f2"o]l that at 
T = 0, C+ remains positive for arbitrarily large fields, 
with (n - 1)C+ « |(1 - x) fo c/& for 6 > & c . 

Thermal effects: Away from 6 C and the XX.Z limit, 
the main thermal effect in Eqs. (|2"2")) - (l2n) will arise from 
the exponential term —2e~@ x , which stems in MF+RPA 
from the Hartree contribution to a M and (s z ) (« jr^/v^ 
and \zjv z ; if just this contribution is kept, Eqs. ©-(|6]) 
lead to C± = i(tanh 2 |/3A— 1) w -2e^ A for /3A > 1). It 
represents the effect of the temperature induced decrease 
of the total spin average (S 2 ) ss jn 2 tanh 2 (3X/2. In the 
exact result it arises from the lowest state of the S — 
n/2 — 1 multiplet, which has excitation energy s» A (see 
Fig. 2]in next section) and multiplicity n — 1 in Eq. ©. 

The RPA thermal factor coth \ fiio cannot, however, be 
neglected (i.e., replaced by 1) in (|2"2"j) - (|2l)1) . particularly 
for b close to b c or \ close to 1, as a; is lower than A (for 
Vfj, > 0) and vanishes for b —¥ b c or v y — > v x . In the exact 
result it represents essentially the effect of the excited 
states within the S = n/2 multiplet, whose excitation 
energies have an approximate harmonic behavior (i.e., 
AE w kw, k = 0, 1, . . .; see Fig. gj}. With this factor, 
Eq. (|23l correctly reduces for u a — » v x and up to 0(l/n) 
terms, to the asymptotic result for the XXZ case 23], 



n-V 1 - (6/6 c ) 2 



while for 6 — ?> 6 C , Eqs. (|22|) and (|25[) converge to 



(29) 



1 9T 

C + « -(1 )-2e-^, (6 = 6 C ). (30) 

n - 1 v x -v y 



Hence, in these regions the concurrence will initially ex- 
hibit an almost linear decrease with increasing T before 
the exponential term becomes appreciable, as a conse- 
quence of the low excitation energy of the S = n/2 states. 

In any case, for sufficiently large n, the concurrence 
will decrease monotonously with increasing T, with C± 
vanishing at a limit temperature that will decrease 
logarithmically with increasing n, as implied by Eq. (|27l) : 



A 



In- 



2(ra-l) 



(31) 



)±! coth ~/3u) 



which is actually a transcendental equation for T^. Both 
vanish (logarithmically) for b — > b f , with T£ decreas- 
ing and increasing with increasing field (and T^(b) 
developing a slope discontinuity at b c ). The increase of 
T£ with increasing b persists for b 3> b c , where 



4(n-l)b/b c 



In 



(32) 



i-x 



implying that at any fixed T parallel entanglement can 
be induced by increasing the field (for & > 6 C , A and uj 
become proportional to 6, the system approaching then 
the entangled ground state as b increases). The same 
behavior was observed in the limit temperatures for non- 
zero global negativities in small anisotropic systems [ljj ■ 
At fixed low T, the main thermal effect for < x < 1 is 
thus the appearance of a separable window 6£ < \b\ < 6 J 
(instead of a separable point) where C± = 0, with 



z y/l- (1 - x)[tanh|/3o;(l - 2(n - l)e~^)] 



±2 



(33) 

(valid for T < T+(b c ) for b+ and T < T£(0) for 67). 
Its width increases then with increasing n or T, with 

< 1. For 



6j ^b s [l±2{x~ 1 -l){n-l)e- fjv *} for' 
T > T£"(0) the separable window will extend through 
b= (C± = for |&| < 6+). 



III. COMPARISON WITH EXACT RESULTS IN 
FINITE SYSTEMS 

Typical results for the magnetic behavior of the con- 
currence at finite temperatures are shown in Fig. [1] for 
the XY case (v z — 0) with n — 100 spins and differ- 
ent anisotropics. It is first seen that MF+RPA results 
obtained with the asymptotic expressions (|22p - (|2"5)) are 
very accurate except in the vicinity of the critical field, 
improving as T increases. The full CSPA results further 
improve those of the MF+RPA in the critical region for 
not too low T, being practically undistinguishable from 
the exact ones at the finite temperatures considered. 

The top panel corresponds to the Ising case v y = 0, 
where the concurrence is always parallel. At T = it 
smoothly increases from as b increases, having a max- 
imum near 6 C , while for T > it becomes non-zero only 



t — i — i — | — i — i — i — | — i — i — i — | — i — i — r 




0.5 1.0 1.5 



b/b c 

FIG. 1. (Color online) Scaled concurrence as a function of the 
magnetic field b for n = 100 spins coupled through a full range 
XY interaction with anisotropies \ — v v /v x = (top), 0.5 
(center) and 0.98 (bottom), at different temperatures. Exact 
and asymptotic mean field+RPA (Eqs. (|22l) <I25|0 results are 
depicted, together with those of CSPA (Eq. |9j for T > 
(almost undistinguishable from the exact ones) . C is parallel 
(antiparallel) for b > b s (< b s ), with b s the factorizing field 
(|28|) . The inset depicts the (parallel) concurrence reentry for 
b > b c at x = 0.98. 

above a threshold field b\ (Eq. ((33j) for T < T^(b c )). In 
the central panel (% = 0.5) we may appreciate the vanish- 
ing of the concurrence at the factorizing field b s ps 0.71u x 
at T = 0, where it changes from antiparallel to parallel. 
This point evolves into a separable window as T increases, 
which extends through b = for T > T£(0) ps 0.15v x . 

The bottom panel depicts the behavior close to the 
XXZ limit. In this case the exact T — concur- 
rence C_ displays an oscillatory behavior as b — > b s 
from below, as in the XXZ chain [23], which reflects 
the ground state spin parity transitions and which is not 
reproduced by MF+RPA (see however discussion of Fig. 
[5]). Nonetheless, as T increases the oscillations become 
rapidly washed out and the asymptotic MF+RPA result 
becomes again accurate, correctly reproducing the exact 
concurrence at T/v x = 0.07 and 0.14, including the reen- 
try of the parallel concurrence that takes place for high 
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FIG. 2. (Color online) Top: Limit temperatures for pairwise 
entanglement Tl as a function of the magnetic field b for 
n — 100 spins at the same anisotropies of Fig. [TJ according to 
exact and MF+RPA results (Eq. IpTC])). They vanish at the 
factorizing field b s . Regions below the limit temperature have 
finite pairwise entanglement, of antiparallel (parallel) type if 
b < b s (> b s ). Bottom: Limit temperatures for increasing 
number of spins at x — 0-5 (n = 10 fc , k = 2, . . . ,6). 



fields. The thermal RPA factor coth i/3w is here essential 
for the accuracy as (3w is small {oj/v x < 0.14 for b < b c ). 

Let us mention that for x € (Qj l]j the ground state, 
which has definite spin parity P = ±1, exhibits n/2 tran- 
sitions ± — > + as b increases from 0, the last one at 
the factorizing field b s . The ground state concurrence 
changes from antiparallel to parallel just at this last tran- 
sition. These transitions are, however, appreciable only 
for x close to 1 (and hence b s close to b c ) or for small sizes, 
as otherwise the ground states of both parity sectors are 
practically degenerate and the concurrence is nearly the 
same in both states (see Fig. 4) as well as in their mix- 
ture. Another consequence of parity conservation is that 
the exact side limits of C± at b = b s are actually non- 
zero and different in finite chains (nC± —> 5/{e s / 2 ± 1), 
with S — n(l — x) [Hj])j being then appreciable for small 
finite (5. In the bottom panel we thus obtain the side 
limits nC- ps 1.16, nC+ ps 0.54 for the exact result at 
b s « 0.996 c , with C± being in fact maximum at b = b s . 

Fig. [2] depicts the corresponding limit temperatures 
T^, which, remarkably, are also accurately reproduced 
by the asymptotic MF+RPA result obtained from Eq. 
(f3"Tj) . Tl vanishes at b = b s but increases V b > b s , devel- 
oping thus a separable field window between the antipar- 
allel and parallel concurrences. The bottom panel shows 
the logarithmic decrease of Tl with increasing n in all re- 
gions. Let us also remark that the behavior of Tl bears 
no relation with that of the mean field critical tempera- 
ture, Eq. (|17[) . which does not depend on the anisotropy 



FIG. 3. (Color online) Thermal behavior of the concurrence 
for x — 0.5 and n = 100 at the indicated fields. The inset 
depicts the non-nonotonous behavior just above b a for \ = 
0.98 and n = 100. Solid lines depict exact results, dotted 
lines those from MF+RPA, Eqs. (p2]) - ([25|l . 



X and vanishes for b > b c . For |&| < b c it is essentially 
higher than Tl (except for very low n [l7j). decreasing 
monotonously from ~v x at b — to at b — b c . 

For x = 0.98, the exact limit temperature T7 actu- 
ally exhibits a small positive slope close to b s , as seen 
in the top panel (not reproduced by MF+RPA). This 
entails that at low finite T the antiparallel concurrence 
will persist in a narrow region above b s , while at fixed b 
within this region, the thermal behavior of the concur- 
rence will be non-monotonous, being first parallel, van- 
ishing and becoming then antiparallel before extinguish- 
ing at the final T7 , as depicted in the inset of Fig. [3] 
Roughly, for small b — b s > 0, it is possible to show that 
T± « cr x (b - b s )5e~ 5 / 2 /{l - e~ 5 ), with a = lncoth<5/4, 
this effect being then noticeable for finite S — n(l — x)- 

The different thermal response of C for fields below, 
around and above the critical field b c can be seen in the 
main panel of Fig.|3]for x = 0.5. The more rapid decrease 
with increasing T for b w b c is in agreement with Eq. (|30[) , 
while the results at b/b c — 0.5 and 1.5 reflect the different 
decrease rate (Eqs. (|23)) - ([25j) ). 

The origin of the distinct thermal factors in Eqs. (122[) 
(f2l)j) can be seen in Fig. 21 which depicts the excitation 
energies of the lowest levels (S, k ± ) for \ = 0.5. For 
\b\ < b c , corresponding levels of opposite parity are prac- 
tically degenerate. The non-vanishing excitation energies 
within the maximum spin multiplet are nearly harmonic, 
the lowest one practically coinciding with the RPA en- 
ergy (fT9|h whereas the excitation energy of the lowest 
state with S = n/2 — 1 is almost b independent and co- 
incident with A = v x . The parity degeneracy becomes 
broken in all levels as b approaches b c , where the max- 
imum spin excitations become low and give rise to the 
increased thermal sensitivity (Eq. (15U|) ). For |6| > b c the 
RPA energy represents again the lowest excitation energy 
in the maximum spin multiplet, which has now negative 
parity, whereas A = b + v z is again the excitation of the 
lowest S = n/2 — 1 state, which has now positive parity. 

The corresponding concurrences are shown in the lower 




FIG. 4. (Color online) Top: Lowest excitation energies 
AE Sk ± = (E Sk ± - E )/v x for X = 0.5 and n = 100, vs. 
magnetic field b. E$k v denotes the energy of level k with to- 
tal spin S and parity v, while Eo = E n / 2 ,a+ ■ Solid (dashed) 
lines depict levels of positive (negative) parity, practically de- 
generate for b < b c . The dotted lines depict the mean field 
and RPA energies A and uj respectively. Bottom: The cor- 
responding concurrences, including that in the ground state. 
For b < b c C is antiparallel (parallel) for b < b s (> b 3 ) in both 
states (n/2, 0*), whereas for b > 6 C , it is parallel in the ground 
state (n/2,0 + ) but antiparallel in the other levels depicted. 



panel. For |6| < b c and x = 0.5, C is non-zero just in the 
degenerate ground states, being almost coincident except 
for b close to b c and changing both from antiparallel to 
parallel at b s . However, for b > b c C is non-zero in all 
maximum spin states, being parallel in the ground and 
highest states but antiparallel in the rest. They are es- 
sentially the basic states IS = n/2, S z = M) plus pertur- 
bative corrections. For \M\ < n/2 they are already en- 
tangled and exhibit hence antiparallel concurrence [23| . 
while for \M\ = n/2 the concurrence arises just from the 
corrections and is hence parallel. We also note that the 
mixture of the n — 1 states with lower spin S = n/2 — 1 
has zero concurrence at all fields (the same occurs with 
lower spin mixtures) so that it can only decrease the ther- 
mal concurrence, which arises then essentially from the 
ground state, except in anomalous regions (the antipar- 
allel reentry in the inset of Fig. [3] arises from the first 
excited state). 

Finally, Fig. [S] depicts results for a small chain (n — 
10), where finite size effects become exceedingly impor- 
tant. The stepwise behavior of the exact concurrence at 
T = is now visible already for x = 0.5 (top panel), the 
side limits at the exact factorizing field bf = (1 — n' r )b s 
being non-zero. The MF+RPA result is now less accurate 
at T = 0, and leads to the vanishing of C_ at a field lower 
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FIG. 5. (Color online) Magnetic behavior of the concurrence 
for a chain with n — 10 spins at different temperatures, for 
anisotropies x = 0.5 (top) and \ — 0.98 (bottom). In the 
latter the antiparallel concurrence increases with increasing 
field at very low T, as in the XX Z case (see text). 

than b s (and closer to b e a x ) if the full expression (IA.2I) is 
used for C_. Nonetheless, MF+RPA results rapidly im- 
prove as T increases, while the CSPA, although no longer 
exact, improves again results in the vicinity of b c . 

When the full square root in the evaluation of C_ is 
kept (Eq. (fA~2jl). the MF+RPA is actually able to qual- 
itatively account for a finite jump in C for fields close to 
b s but only if \ is very close to 1, as seen in the bottom 
panel. In this case, the MF+RPA result for C_ at T = 
does not decrease as b increases but rather increases, in 
agreement with the exact result, terminating at a final 
field bf < b e s x where it starts to be complex (and is max- 
imum). For x = I — S/n and fields just below b c , i.e., 
(b/b c ) 2 = 1 — e/n, we actually obtain, instead of Eq. 
([23]). the asymptotic MF+RPA expression 

C_ » i(I £ - y^coth \(3lu) - 2e-0"- 

-ip^/fcothi^ + iECe-^] 1 ^ (34) 

with uj = \fe5(v x — v z )/n. While for e oc n it reduces to 
Eq. ([231) j for e < 4 it becomes complex if 6 is sufficiently 
small. At T = 0, if S < S c = ^ 0.55, Eq. (|3i|) 
becomes complex for e < £f(5) ps 2.4 + | V<5 C — (5, with 
C_(6y) « Note that c_(b f ) > 1 for 6 < 0.48, with 
C-(bf) — > 2/n for 5 — > 0, which is the correct result for 
the 112 limit [H). In the case depicted, 5 = 0.2. 

Eq. implies as well that these effects, in particular 
the increase of C_ as b — > b s for \ close to 1, will dis- 
appear for very low T oc uj oc 6 c /n, which is confirmed 



in the exact results. One may also appreciate in the 
bottom panel of Fig. [5] the significant persistence of the 
antiparallel concurrence up to b w 1.56 c at T/v x = 0.14 
(not reproduced by MF+RPA or CSPA), which is just 
the same anomalous behavior discussed in Figs. [2]-[3l en- 
hanced by the smaller value of S. Nonetheless, even in 
this extreme case there is a weak but non-zero revival of 
the parallel concurrence for high fields V T (appreciable 
in the figure just for T = 0) which is correctly reproduced 
by MF+RPA. 



IV. CONCLUSIONS 

We have analyzed the thermal behavior of the pair- 
wise concurrence in a fully connected spin system with 
anisotropic couplings placed in a transverse field. For 
the usual 1/n scaling of coupling strengths, the limit 
temperature decreases only logarithmically as the size 
n increases, and decreases (increases) for increasing field 
in the antiparallel (parallel) sectors, the latter extend- 
ing for arbitrarily large fields. This behavior was previ- 
ously observed for small n in the temperatures limiting 
global negativities, for which the pairwise limit tempera- 
ture provides a lower bound. Anisotropic arrays become 
then strictly pairwise separable just within a finite field 
window at any temperature, which collapses into the fac- 
torizing field at T = 0. Remarkably, all previous features 
of the pairwise entanglement can be captured by a simple 
thermal MF+RPA treatment, consistently derived from 
the path integral representation of the partition function, 
which in the present case is able to provide a reliable an- 
alytic description of the concurrence and limit tempera- 
ture at all fields, exact in the large n limit. We have also 
discussed the special finite size effects exhibited by this 
system for small anisotropies or sizes, whose main as- 
pects can also be qualitatively reproduced by MF+RPA 
or the full CSPA. These results suggest the possibility 
of describing by simple means at least the main features 
of the thermal pairwise entanglement in more complex 
systems, although the actual accuracy and scope of the 
RPA in such situations remains to be investigated. 

The authors acknowledge support from CIC(RR) and 
CONICET (JMM, NC) of Argentina. 



Appendix: MF+RPA Concurrence 

We provide here the expressions for the full MF+RPA 
concurrence derived from Eqs. ©-([S]) and ([i"6]) . Setting 
v x = 1 and b = b/b c , in the symmetry breaking phase we 
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obtain For /3A < 1, A w 1 - 2e~ /3A , with 1 - A 2 « 4e _,3A . Eqs. 

22 ]> -([23 ]l are then obtained from (jA.ljl . (|A.3|) neglecting 



C+ = 1 {1 — — cothi/3w[l + C an d setting A = 1 in the 0(l/n) terms. This is cor- 

^ n ~ ^ ~ v y rect up to 0(1/ n) terms for temperatures where C± are 

£ (1 — v y )\ 2 ^ ( 2 . C)T]} ^ ^ positive, as in such a case e~^ A must be 0(l/n). 
1 — C A 2 — b 2 1 — C 
A 2 -6 2 1 r l-6 2 

1 -% coth i^ [ ^+^+ c a2(1 _ v 



Similarly, in the normal phase we obtain 



^ u "2^~l 2 1-C" V " " zn l-tanh 2 i/3A 1 r , 
c2 C + = ^ + _ { i__^ ooth ^ w + 



(i-C) 



[l-(3-C)T]} 



11 2 n-l v w 2 2 2 n ( 3T 

_^ (1 + i^ coth i^ )2} i/ 2) (A2) ~(w)~ } |A ' 4) 

nw 



where £ = |/3(1 — A 2 ) and A = tanh If b is not close 
to 1, we may expand Eq. (|A.2|) up to 0(1/ n) as 

^ 1-A 2 1 . , . r A 2 -& 2 and A determined by Eq. (EOl). For 8X < 1, 1 - 

^ + :-^U-— ^cothi/? w [— ^- . ,21, 



where / M = v M tanh(±/3A)/A, C = \Pv z (l - tanh 2 |/3A) 



2 + n - l 1 w 2P [ l - &2 tanh 2 ±/3A w -4 e -' 3A . Eq. (Q leads then to Eq. (2£ 

a ^2 up to 0(l/n) for temperatures where C+ > 0, by ne- 

(1 - u 2 )A 2 ] - - [1 - (3 - C)r](}A.3) glccting \ and setting \ — b + v z in the 0(1/ n) terms. 
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